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The probability density function of single-point velocity fluctuations in turbulence is studied 
systematically using Fourier coefficients in the energy-containing range. In ideal turbulence where 
energy-containing motions are random and independent, the Fourier coefficients tend to Gaussian 
and independent of each other. Velocity fluctuations accordingly tend to Gaussian. However, if 
energy-containing motions are intermittent or contaminated with bounded-amplitude motions such 
as wavy wakes, the Fourier coefficients tend to non-Gaussian and dependent of each other. Velocity 
fluctuations accordingly tend to non-Gaussian. These situations are found in our experiment of a 
rough-wall boundary layer. 

PACS numbers: 47.27.Ak, 47.27.Nz 



I. INTRODUCTION 

Suppose that single-point velocity fluctuations u(x) are 
measured repeatedly in stationary turbulence over the 
range < x < L. If the sampling interval is much greater 
than the eddy turnover timescale, these measurements 
serve as independent realizations of the turbulence. Each 
of them is expanded into a Fourier series as 




where 2nn/L — k n is the wave number. Batchelor 
assumed that the Fourier coefficients a n and b n are sta- 
tistically independent of each other and applied the cen- 
tral limit theorem to their sum. This theorem ensures 
that the probability density function (PDF) of a sum of 
many independent random variables tends to Gaussian, 
at least within a few standard deviations around the av- 
erage 0, ■ It was concluded that velocity fluctuations 
tend to Gaussian, being consistent with experimental and 
observational data of turbulence that were available at 
that time. 

However, recent measurements revealed the presence of 
velocity fluctuations that tend to non-Gaussian. Sreeni- 
vasan and Dhruva Q obtained long data of atmospheric 
turbulence at 35 m above the ground. Their data yield 
the flatness factor F u = (u A ) / (u 2 ) 2 = 2.66, where the 
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bracket (■) denotes an average. This value is different 
from the Gaussian value F u — 3. Noullez et al. ||| ob- 
tained F u ~ 2.85 in turbulent jets. For these results, 
there has been no explanation. 

Therefore, although single-point velocity fluctuations 
are fundamental in describing turbulence, the mecha- 
nism that determines their PDF is uncertain. We sys- 
tematically study velocity fluctuations in a laboratory 
rough-wall boundary layer, a representative turbulent 
flow with various applications, e.g., the atmosphere near 
the ground. With an increase of the distance from the 
wall, velocity fluctuations are found to change from sub- 
Gaussian (F u < 3) to Gaussian, and to hyper-Gaussian 
(F u > 3). This behavior is discussed using PDFs of the 
Fourier coefficients and correlations among them. 



II. CONDITION FOR GSAUSSIANITY 

This section serves as a summary of conditions for ve- 
locity fluctuations to be approximately Gaussian. It is 
assumed that the turbulence is not only stationary but 
also homogeneous in the x direction. The data length 
is set to be much larger than the correlation length 
l c = j \{u(x + 5x)u(x))\dSx I (u 2 ) . Then an average taken 
over the realizations is equal to the corresponding aver- 
age taken over the x positions. 

Velocity fluctuations of turbulence are dominated by 
Fourier coefficients in the energy-containing range. Since 
the data length is large, there is a sufficient number of 
Fourier coefficients for the central limit theorem to be 
applicable. Velocity fluctuations tend to Gaussian if the 
Fourier coefficients are random and independent. This 
is expected for ideal turbulence where energy-containing 
motions are random and independent. Although tur- 
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bulence always contains small-scale coherent structures, 
e.g., vortex tubes their contribution to velocity fluc- 
tuations is as small as the energy ratio of the dissipation 
range to the energy-containing range. 

The central limit theorem is not applicable to a sum of 
variables if few of them dominate over the others 0, 01 • 
For example, if the energy spectrum E n = (a^) + (b^) 
is proportional to k" with a < — 1, Fourier coefficients 
at the smallest wave numbers dominate the velocity fluc- 
tuations. They do not necessarily tend to Gaussian 0. 
Nevertheless, the energy spectrum is relatively flat in the 
usual energy-containing range, where the power law fc" 
is not a good approximation. We previously assumed 
that the central limit theorem is not applicable to the 
sum of Fourier coefficients in any turbulence , but this 
assumption was wrong. 

Fourier coefficients of velocity fluctuations also tend 
to Gaussian if turbulence is made of random and inde- 
pendent motions [aSllfl- The Fourier coefficient a n is 
obtained as 
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For n 1, we are able to set l<m<n and l c <C L/m. 



The integrations ... dx, and f/' ,x r/ . 

& JO ' ' J (m — l)L/m 

regarded as random variables. Their amplitudes are the 
same because the segment size L/m exceeds the wave- 
length L/n. They are nearly independent because adja- 
cent integrations depend on each other only at the ends. 
Then a n tends to Gaussian as a consequence of the cen- 
tral limit theorem. This discussion is not applicable to 
Fourier coefficients a n and b n for n ~ 1. They are nev- 
ertheless small and do not contribute to velocity fluctua- 
tions if the data length L is sufficiently large, because of 
the universal trend E n — > in the limit k n — 2im/L — > 
0. We are able to assume safely that all the Fourier 
coefficients are Gaussian. 

The central limit theorem offers no information about 
the tails of the PDF of a sum of variables. Only when 
the sum (minus its mean) has been divided by the square 
root of the number of the variables, the Gaussian approx- 
imation holds at the tails 0, HI • Although statistics such 
as the flatness factor mainly reflect the core of the PDF, 
the tails that are significantly far from Gaussian could 
exist and affect the statistics. One example is significant 
contamination with a bounded-amplitude motion, e.g., 
u{x) oc sin(x), the PDF of which does not have tails. Ve- 
locity fluctuations tend to sub-Gaussian, e.g., F u = 3/2 
for u{x) oc sin(a;) 0. The Fourier coefficients that corre- 
spond to the bounded-amplitude motion also tend to sub- 
Gaussian [Tlj . Another example is significant intermit- 
tency, where velocity fluctuations tend to hyper-Gaussian 
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FIG. 1: Schematic representation of our experimental setup. 
Plan and side views are shown together with coordinate axes 
(upper and lower sketches, respectively). 



.12]. The Fourier coefficients also tend to hyper-Gaussian 
because some of the integrations J^ m ...dx, and 
f, L T i ...dx are enhanced over the others. Since the 

J \m—\)L/m 

Fourier coefficients are linear functions of velocity fluctu- 
ations, sub- and hyper-Gaussianities of the velocity fluc- 
tuations yield sub- and hyper-Gaussianities of the Fourier 
coefficients, respectively, and vice versa. We expect that 
this relation could approximately exist even in the case 
of mutually dependent Fourier coefficients as long as the 
dependence is weak. 



III. EXPERIMENT 

The experiment was done in a wind tunnel of the Mete- 
orological Research Institute. As shown in Fig. 1, we use 
the coordinates x, y, and z in the streamwise, spanwise, 
and floor- normal directions. The corresponding wind ve- 
locities are u, v, and w. We take the origin x = y = z = 
on the tunnel floor at the entrance to the test section. 
The test section had the size 5x = 18 m, Sy — 3 m, and 
5z = 2 m. A boundary layer was made by placing blocks 
over the entire floor of the test section. The blocks had 
the size dx = 6 cm, Sy = 21 cm, and 5z = 11 cm. The 
spacings of the adjacent blocks were 5x = 70 cm and 
5y = 80 cm. We set the in-coming wind velocity to be 20 
m s _1 . The boundary layer was well developed at x > 10 
m. 

The u and v or u and w components of the wind ve- 
locity were measured using a hot-wire anemometer with 
an X-type probe. The wires were made of tungsten, 5 
/im in diameter, 1.0 mm in effective length, 1.4 mm in 
separation, oriented at ±45° to the streamwise direction, 
and operated at the temperature of 280°C. The measure- 
ment positions were at x = 10 m and z — 0.05-1.00 m. 
The signal was low-pass filtered with 24 dB/octave and 
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FIG. 2: (a) Mean streamwise velocity U. (b) Root-mean- 
square velocity fluctuations (u 2 ) 1 / 2 , (v 2 ) 1 / 2 , and (w 2 ) 1 ^ 2 . (c) 
Shear-stress velocity (—uw) 1 ^ 2 . The abscissa is the height z. 
The triangles denote the u component, the circles denote the 
v component, and the squares denote the w component. The 
horizontal dotted lines indicate the values of U and (v 2 ) 1 ^ 2 ob- 
tained in the in-coming flow. The arrows indicate the height 
of the roughness elements and the displacement thickness. 
On the uppermost axis, we indicate the height normalized by 
the 99% thickness. On the right axis, we indicate the val- 
ues normalized by the maximum U value. We also indicate 
the z range of the constant-stress sublayer. The shear-stress 
velocity is not available at z — 0.90 and 1.00 m, where the 
correlation —uw is negative. 



sampled digitally with 16-bit resolution. At z < 0.30 m, 
the filtering frequency was 10 kHz and the sampling fre- 
quency was 20 kHz. At z > 0.30 m, they were 25 kHz 
and 50 kHz, respectively. The length of the signal was 
8 x 10 6 points at z = 0.11, 0.14, 0.24, 0.28, 0.40, 0.50, 
0.60, 0.80, and 1.00 m or 32 x 10 6 points at z = 0.05, 
0.08, 0.17, 0.20, 0.32, 0.36, 0.70, and 0.90 m. We used 
the frozen-eddy hypothesis of Taylor to convert temporal 
variations into spatial variations. 

Figure [21 shows the mean streamwise velocity U, the 
root-mean-square velocity fluctuations (u 2 ) 1 ^ 2 , (v 2 ) 1 ^ 2 , 
and (w 2 ) 1 / 2 , and the shear-stress velocity (—uw) 1 ' 2 as 
a function of the height z. The 99% thickness, i.e., the 
height at which U is 99% of its maximum value U, is 0.68 
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FIG. 3: (a) Flatness factors F u , F v , and F w . (b) Skewness 
factors S u , S v , and S w . The abscissa is the height z. The 
triangles denote the u component, the circles denote the v 
component, and the squares denote the w component. The 
horizontal dotted lines indicate the Gaussian values F = 3 
(a) and S = (b). On the uppermost axis, we indicate the 
height normalized by the 99% thickness. We also indicate the 
z range of the constant-stress sublayer. 



m. The displacement thickness f Q (1 — U/U)dz is 0.17 
m For reference, the average U and the root-mean- 
square fluctuation (i* 2 ) 1 / 2 obtained in the in-coming flow 
at x — z — 1 m are shown. The latter reflects the me- 
chanical and electric noise, which is well below the tur- 
bulence signals. 

Throughout the boundary layer, velocity fluctuations 
are almost isotropic [Fig.f^b)]. Although velocity fluctu- 
ations in a smooth-wall boundary layer are anisotropic, 
the anisotropy is reduced over roughness [l4| . The shear- 
stress velocity is almost constant at z ~ 0.14-0.40 m 
[Fig. life)]. Below and above this constant-stress sub- 
layer, there are the roughness sublayer and the outer 
sublayer, respectively, where turbulence is affected by the 
roughness or the outer laminar flow. The logarithmic law 
that corresponds to the constant stress is unclear in the 
profile of the mean streamwise velocity [Fig. GJa)] be- 
cause our constant-stress sublayer is relatively thin and 
there is uncertainty in defining the effective origin for the 
logarithmic law. 
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IV. RESULTS AND DISCUSSION 

Figure Ofa) shows the flatness factors F u , F v , and F w 
as a function of the height z. Figure OJb) shows the 
skewness factors S u — {u 3 } / (u 2 ) 3 ^ 2 , S v , and S w . The 
flatness factors are close to 3 at z < 0.40 m. As the height 
is increased above z ~ 0.40 m, the flatness factors begin 
to increase. They have pronounced peaks at z ~ 0.90 m. 
The skewness factors are also significant at z > 0.40 m, 
except for the v component that is free from the shear 
of the boundary layer. Similar results were obtained in 
previous works [TH Hil . 

We focus on the v component, which is best suited to 
our analysis. The range of flatness factor is widest. The 
dependence on the height z is simple. At z < 0.40 m, 
the PDF is sub-Gaussian. At z > 0.40 m, the PDF is 
hyper-Gaussian. Only the v component at z ~ 0.40 m 
exhibits simultaneously the Gaussian values F = 3 and 
5 = 0. 

Figure Ufa) shows the energy spectrum E n for the 
v component at z = 0.05, 0.36, and 0.70 m. Fig- 
urcUb) shows the flatness factor of the Fourier coefficient 
F n = (a*)/(a 2 ) 2 . FigureHfc) shows the mode-mode cor- 
relation C n n p (n ^ Tip): 
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where n p corresponds to the wave number k n of the 
E n peak. These quantities were obtained by dividing the 
data into segments of 2 15 points. We regarded them as in- 
dependent realizations of turbulence, applied the Fourier 
transformation individually to them, and calculated aver- 
ages over them at each of the wave numbers. The energy 
spectrum was obtained using the Welch window function. 
Since this and other usual window functions were found 
to affect significantly the flatness factor or mode-mode 
correlation, we were forced to obtain them by appending 
the inverted sequence to each sequence of the segments. 
Any method to remove effects of discontinuity at the ends 
of a data segment modifies the data and thereby affects 
some of the statistics. Our present method is not an 
exception but happens to have no serious effect on the 
statistics of our interest. 

Velocity fluctuations at z = 0.36 m are Gaussian 
[F v = 3.00; Fig. Efa)]. Throughout the wave numbers, 
the flatness factor of the Fourier coefficient is close to 
Gaussian [Fig. 0Jb)], and the mode-mode correlation is 
absent [Fig.^Jc)]. Thus the Fourier coefficients are Gaus- 
sian and independent of each other. The height z = 0.36 
m is near the upper edge of the constant-stress sublayer 
and also near the middle of the entire boundary layer. 
We consider that eddies of various sizes and strengths 
pass the probe randomly and independently. 

Velocity fluctuations at z = 0.05 m are sub-Gaussian 
[F v = 2.51; Fig. Qta)]. At around the peak of the energy 
spectrum E n , the flatness factor of the Fourier coeffi- 
cient is sub-Gaussian [Fig. Bfb)]. This sub-Gaussianity 
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FIG. 4: (a) Energy spectrum E n of the v component at 
z = 0.05, 0.36, and 0.70 m. (b) Flatness factor of the Fourier 
coefficient F n . (c) Mode-mode correlation C nnp (n ^ n p ). 
The abscissa is the wave number k n . The horizontal dotted 
lines indicate the Gaussian value of 3 (b) and the noncorrela- 
tion value of (c). The arrows indicate the wave numbers k np 
of the En peaks. In (b) and (c), we made moving averages 
over five adjacent wave numbers. The hatched areas empha- 
size the sub-Gaussianity and negative correlation of the data 
at z = 0.05 m. 



of the Fourier coefficients is associated with that of ve- 
locity fluctuations (Sec.[nJ). Although the former is less 
significant than the latter, the mode-mode correlation is 
negative [Fig.^c)]. This is also associated with the sub- 
Gaussianity of velocity fluctuations. Even if the Fourier 
coefficient at the E n peak has a large amplitude, its effect 
to velocity fluctuations tends to be weakened by small 
amplitudes of the Fourier coefficients at nearby wave 
numbers, as compared with a noncorrelation case. 

The mode-mode correlation shown in Fig. E^c) is 
merely a representative example. Similar negative cor- 
relations are found for other pairs of Fourier coefficients 
at around the energy peak, ft is possible to have more 
than two Fourier coefficients that are negatively corre- 
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lated with each other if the absolute values of the corre- 
lation coefficients are small. 

The height z = 0.05 m is in the roughness sublayer. A 
plausible explanation is that turbulence is contaminated 
with bounded-amplitude motions due to wavy wakes of 
the roughness. The amplitudes are required to be in a 
bounded range because velocity fluctuations and their 
Fourier coefficients are sub-Gaussian (Sec. ITT)!. The indi- 
vidual motions are required to contribute to a range of 
wave numbers, possibly through the presence of spatial 
structures, because there are mode-mode correlations. 

Velocity fluctuations at z = 0.70 m are hyper-Gaussian 
[F v = 14.82; Fig. Ola)]. Throughout the wave num- 
bers, the flatness factor of the Fourier coefficient is hyper- 
Gaussian [Fig. Ub)]. This is associated with the hyper- 
Gaussianity of velocity fluctuations (Sec.[nJ). The mode- 
mode correlation is positive [Fig.^c)]. If the Fourier co- 
efficient at the E n peak has a large amplitude, its effect to 
velocity fluctuations is strengthened by large amplitudes 
of the Fourier coefficients at nearby wave numbers. 

The height z = 0.70 m is near the outer edge of the 
boundary layer, where turbulence is intermittent [l2lll7j|. 
There are only eddies that have been ejected from the 
lower heights. Actually, the skewness factors of the u 
and w components are negative and positive, respectively 
[Fig. GUb)]. The eddies intermittently pass the probe 
and enhance velocity fluctuations. They are accordingly 
hyper-Gaussian. The intcrmittcncy is also responsible for 
the hyper-Gaussianity of the Fourier coefficients fSec. lTTjl . 
The presence of spatial structures is responsible for the 
presence of mode-mode correlation. It should be noted 
that, since turbulence is not space filling, this intermit- 
tency is regarded as contamination of turbulence with a 
laminar flow. 

The hyper-Gaussianity of velocity fluctuations in the 
outer sublayer has been known to be an intermittent phe- 
nomenon for a long time, readily understood in the space 
domain Their Fourier coefficients in the wave num- 
ber domain are not so useful but have interesting proper- 
ties as in the cases of sub-Gaussian and Gaussian velocity 
fluctuations. 

Figure El summarizes behaviors of the Fourier coeffi- 
cients as a function of the height z: (a) the wave number 
k np of the peak of the energy spectrum, (b) the flatness 
factor of the Fourier coefficient F n , and (c) the mode- 
mode correlation C nn (n ^ n p ). The latter two quan- 
tities are medians around k 7lp within 2 ±1 fc Ilp . Since sta- 
tistical uncertainty is not insignificant, the shorter data 
with 8 x 10 6 points are denoted by the smaller symbol. 
With an increase of the height, the peak wave number de- 
creases because energy-containing eddies become larger 
[l8|. The behaviors of the flatness factor and mode- 
mode correlation are in accordance with the behavior of 
velocity fluctuations. At z < 0.40 m, velocity fluctua- 
tions are sub-Gaussian. The Fourier coefficients are sub- 
Gaussian and exhibit negative mode-mode correlations. 
At z > 0.40 m, velocity fluctuations are hyper-Gaussian. 
The Fourier coefficients are hyper-Gaussian and exhibit 
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FIG. 5: (a) Wave number k np of the peak of the energy spec- 
trum of the v component, (b) Flatness factor of the Fourier 
coefficient F n . (c) Mode-mode correlation C nnp (n 7^ n p ). 
The flatness factor and the mode-mode correlation are me- 
dians around k„ p within 2 ±1 k rlp . The abscissa is the height 
z. The large circles denote data with 32 x 10 6 points, while 
the small circles denote data with 8 x 10 6 points. The hori- 
zontal dotted lines indicate the Gaussian value of 3 (b) and 
the noncorrelation value of (c). On the uppermost axis, 
we indicate the height normalized by the 99% thickness. We 
also indicate the z range of the constant-stress sublayer as 
well as the z ranges where velocity fluctuations are sub- and 
hyper-Gaussian. 



positive mode-mode correlations. 



V. CONCLUDING REMARKS 

The PDF of single-point velocity fluctuations in tur- 
bulence is not universal but reflects energy-containing 
motions. We studied velocity fluctuations using their 
Fourier coefficients in the energy-containing range. In 
ideal turbulence where energy-containing motions are 
random and independent, the Fourier coefficients tend to 
Gaussian and independent of each other. Velocity fluctu- 
ations tend to Gaussian. This is the case at around the 
middle of a rough-wall boundary layer, where eddies of 
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various sizes and strengths pass the probe randomly and 
independently. However, if turbulence is contaminated 
with bounded-amplitude motions such as wavy wakes, 
the Fourier coefficients tend to sub-Gaussian and their 
amplitudes are correlated negatively. Velocity fluctua- 
tions tend to sub-Gaussian. This is the case in the lower 
part of a rough-wall boundary layer, where contamina- 
tion with wavy wakes of the roughness is significant. If 
turbulence is intermittent or contaminated with a lami- 
nar flow, the Fourier coefficients tend to hyper-Gaussian 
and their amplitudes are correlated positively. Velocity 
fluctuations tend to hyper-Gaussian. This is the case in 
the upper part of a rough-wall boundary layer, where 
turbulence is not space filling and eddies intermittently 
pass the probe. 

We previously studied velocity fluctuations in grid tur- 
bulence 8] . At small distances from the grid, turbulence 
is developing. There are bounded-amplitude motions due 
to wavy wakes of the grid. The PDF of velocity fluctua- 
tions tends to sub-Gaussian. At intermediate distances, 
turbulence is fully developed. The PDF tends to Gaus- 
sian. At large distances, turbulence is decaying. There 
remain only strong eddies. They intermittently pass the 
probe. The PDF tends to hyper-Gaussian. These results 
for grid turbulence are consistent with our present results 
for boundary-layer turbulence. 

Sreenivasan and Dhruva Q obtained F u — 2.66 in the 
atmospheric boundary layer. The exact observational 



condition is unknown to us, but the flow at the measure- 
ment position z = 35 m could be affected by wavy wakes 
of the ground roughness. The observed sub-Gaussianity 
could be attributable to possibly bounded amplitudes of 
these wakes. 

Velocity fluctuations tend to sub-Gaussian in direct 
numerical simulations of homogeneous, isotropic, station- 
ary turbulence 0, |2(], . Since the simulations were 
done under forcing over narrow ranges of the smallest 
wave numbers, the energy spectra are steep and close 
to the power law with a < — 1. The observed sub- 
Gaussianity could be attributable to the forced motions 
that dominate the velocity fluctuations Q, El ( see a l so 
Sec.[nJ. Thus these numerical results are not inconsistent 
with our experimental results where the energy spectrum 
is relatively flat in the energy-containing range [Fig.Qfa)]. 
It is of interest to study such numerical data in the same 
manner as in the present work. 
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